Universal behavior of the Shannon mutual information of critical quantum chains 
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We consider the Shannon mutual information of subsystems of critical quantum chains in their 
ground states. Although this quantity depends on the basis the ground state is expressed our results 
indicate that its leading behavior is basis independent. Moreover, as happens with the entanglement 
entropy, its finite-size behavior yields the conformal anomaly c of the underlying conformal field 
theory governing the long distance physics of the quantum chain. We studied analytically a chain of 
coupled harmonic oscillators and numerically the Q-state Potts models (Q — 2, 3 and 4), the XXZ 
quantum chain and the spin-1 Fateev-Zamolodchikov model. The Shannon mutual information is a 
quantity easily computed, and our results indicate that for relatively small lattice sizes its finite-size 
behavior already detects the universality class of quantum critical behavior. 
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Entanglement measures have emerged nowadays as 
powerful tools for the study of quantum many body 
systems [l], In one dimension, where most quantum 
critical systems have their long-distance physics ruled by 
a conformal field theory (CFT), the entanglement en- 
tropy has been proved the most important measure of 
entanglement. It allows one to identify the distinct uni- 
versality classes of critical behaviors. Let us consider a 
periodic quantum chain with L sites, and partition the 
system into subsystems A and B of length I and L — I, 
respectively. The entanglement entropy is defined as the 
von Neumann entropy of the reduced density matrix pj^ 
of the partition A: Si = —Trj^pj^lnp^. If the system is 
critical and in the ground state, in the regime where the 
subsystems are large compared with the lattice spacing, 
Si is given by d, [ 



Is, 



(1) 



where c is the central charge of the underlying CFT 
and 7 S is a non-universal constant. A remarkable fact 
is that even in the case where the system is in a pure 
state formed by an excited state, the conformal anomaly 
dictates the overall behavior of the entanglement, simi- 
larly as in ([T]) Q. It is worth mentioning that recently 
many interesting methods were proposed to cal- 

culate the entanglement entropy and ultimately central 
charge, however, up to know they have not been imple- 
mented experimentally. A natural question concerns the 
possible existence of other measures of shared informa- 
tion that, similarly as the entanglement entropy, are also 
able to detect, the several universality classes of critical 
behavior of quantum critical chains. 

In this Letter we present results that indicate that the 
Shannon mutual information is such a measure. The 
Shannon mutual information of the subsystems A and 
B, of sizes I and L — I is defined as 



where Sh(X) = — ^2 x p x ^npx is the Shannon entropy 
of the subsystem X with probabilities p x of being in a 
configuration x. These probabilities, in the case where 
A is a subsystem of a quantum chain with wavefunc- 
tion \^aub) = En,m c ™,™l<%) ® \4>b), are given by the 
square of the normalized amplitudes of the reduced wave- 
function |*^) = E„ (Em <V™ )!<%)> wh ere {|<%)} and 
{\4>B) m } are the vector basis in subspaces A and B. We 
can think about Shannon mutual information as the von 
Neumann mutual information of the mixed state pro- 
duced after local projective measurements @, @. In [l(| 
it was conjectured that the mutual information, like the 
entanglement entropy should follow the area law; see 
also [Tl(. Many authors studied the Shannon entropy 
of the one dimensional quantum spin chains (T2l - fl4j and 
found useful applications in classifying one and two di- 
mensional quantum critical points. Several authors also 
studied different properties of the Shannon mutual in- 
formation in two dimensional classical systems [HI, [l6| . 
They found that although the mutual information of two 
halves of a cylinder has its maximum value at a temper- 
ature higher than the critical temperature, its derivative 
diverges at the critical temperature. Most recently Um 
et al [9( studied the mutual information of a subregion 
with respect to the rest in the periodic transverse Ising 
model and surprisingly found that it has the same depen- 
dence ln(sin(7r£/L)) as ([TJ) but with a distinct multiplica- 
tive constant. Although the mutual information just re- 
cently has found many applications in spin chain studies, 
it is worth mentioning that the concept has been already 
briefly studied in [l7l [l8j in the context of harmonic oscil- 
lators. In this Letter we study the Shannon mutual infor- 
mation in different critical spin chains and argue about 
their possible connections to the central charge of the 
underlying CFT. Our results based on the study of sev- 
eral quantum chains suggest that for periodic chains in 
the ground state, the Shannon mutual information in the 
scaling regime {t, L >> 1) is universal and has a depen- 
dence with the subsystem size £ similar as the one of the 



I(A, B) = Sh{A) + Sh{B) - Sh(A U B), 



(2) 



2 



entanglement entropy, i. c 
c 



I{l,L) 



In 



-sm(— )) + 7J , (3) 



where c is the central charge of the underlying CFT and 
7/ is a non- universal constant. Up to the moment, in con- 
trast with the entanglement entropy, there is no simple 
general field theoretical method to calculate the Shan- 
non entropy and consequently the Shannon mutual in- 
formation. The difficulty comes from the evaluation of 
the summation over the amplitudes of the ground state 
cigenfunction. 

In order to justify our conjecture §3§ we calculate the 
Shannon mutual information for quantum systems in fi- 
nite geometries. We first present our analytical results for 
a system of coupled harmonic oscillators (Klein-Gordon 
field theory) and then our numerical analysis for sev- 
eral critical quantum spin chains: Q-state Potts model 
(Q = 2, 3, and 4), spin-i XXZ spin chain and the spin- 
1 Fateev-Zamolodchikov model in the antiferromagnetic 
and ferromagnetic regime. 

Harmonic oscillator(Klein-Gordon field the- 
ory) - The Hamiltonian of L coupled harmonic oscil- 
lators with coordinates $i, . . . , <&l and conjugated mo- 
menta 7Ti , . . . , 7T£ is given by 



L 

E ■ 

n.n' — l 



(4) 



where in the case of nearest-neighbor couplings the in- 
teraction K matrix is just the discrete Laplacian. In the 
continuum limit the above Hamiltonian is the one of a 
simple scalar free field theory (central charge c = 1). 
Let us now consider $>a = {<f>\, 4>2, ■ ■ ■ , 4>i) and <&b = 
{4>e+ij 4>e+2, ■ ■ ■ , 0l) as the position vectors of the subsys- 
tems A and B and Ha g the respective momentum vec- 
tors. The Shannon mutual information I (A, B) =I(£, L) 
between two regions A and B is 



I(£,L) 



d L $p($ A , $ B ) In 



where £>($4,$e) = |*o| 2 



is the total 
2 and p 2 ($ 



(5) 



and 



I [n m e(A)#m| l*o| 2 are the 
densities in position space (^o(<Pi 
state wave function). Then after simple integrations one 
can get [13] 



reduced probability 
. . , 4>l) is the ground 




FIG. 1: (Color online) 1(1, L) - 
subsystem size for the Q — 2, 3 
the Ri- and Si-basis. 
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FIG. 2: (Color online) I(£,L) - /(■§,£) versus ln(|sin(^)) 
for the Q = 2, 3 and 4 Potts models with different lattice sizes 
L. The dashed straight lines are obtained from the fitting 
by using the largest lattice size shown for each model. (The 
ground states were taken in the i?i-basis). 



(f>i4>j >= \{K and (PX)ij ■=< W^j >= 
We noticed that the above formula is exactly equal to 
the quantum Renyi entanglement entropy with n = 2 
[lH and consequently using the CFT techniques Q one 
can get the following result for a periodic system 



1 T ttP 
^,i) = -ln(-sin(-)), 

4 7T L 



(7) 



I(£,L) = ^lm>« 



(6) 



where vi are the eigenvalues of the matrix C = s/XaPa 
and Xa and Pa are I x £ matrices describing correlations 
of position and momentum within subsystem A fl9| . In 
other words for < i,j < £ + 1 we have (XA)%j '■=< 



that agrees with the conjecture ([3]). It is worth mention- 
ing that the mutual Shannon information I (A, B) ob- 
tained in momentum basis also follows the same formula 
which indicate that the coefficient 1/4 in (J7J is indepen- 
dent of the basis we choose. 

Quantum Q-state Potts model - The model in a 
periodic lattice is defined by the Hamiltonian (2p| 



3 



L Q-l 

h q = -EE (Wi* + x ^ ( 8 ) 

i=l fe=0 

where 5, and i?i are QxQ matrices satisfying the follow- 
ing Z(Q) algebra: [i? l5 i?,] = [S t ,Sj] = [Si,Rj] = for 
i ^ j and SjRj = t ' o l< s S, and R? = S? = 1. The sys- 
tem is critical at the self dual point A = 1. For Q = 2, 3 
and 4 its critical behavior is governed by a CFT with 

central charge c = 1 , 6 , where JQ = 2 cos ^r-r. 

We first calculate the ground state of the Hamiltonian 
© for Q = 2,3 and 4 in different basis by exact di- 
agonalization. Due to the Perron-Frobcnius theorem the 
coefficients of the ground state eigcnfunction are real and 
their square give us the probabilities for the calculation of 
the Shannon entropy Differently from the von Neumann 
entanglement entropy the Shannon entropy depends on 
the particular basis used for computing the probabili- 
ties. We verified for the critical chains we studied that 
as the lattice size increases the Shannon mutual infor- 
mation exhibits a dominant basis-independent behavior, 
apart from an additive basis dependent constant, as con- 
jectured in ([3]). In order to illustrate this result we show 
in Fig. 1 the difference L) — I(L/2, L) obtained from 
the ground-state at the critical point of the Q = 2, 3 and 
4 Potts models. The calculations were done by expressing 
the ground state in the basis where either the matrices Ri 
or Si arc diagonal. Apart from the initial point £ = 1 we 
already see for these lattice sizes a quite good agreement 
among the results in both basis. In order to show the 
dominant ^-dependence of I(£,L) and test © we con- 
sider the difference I(£, L) — I(L/2, L), since in this case 
the basis dependent constant 7/ is canceled. In Fig. 2 
we plot this difference, as a function of \n(sm(ir£/ L)/4, 
for some lattice sizes of the Q — 2, 3 and 4 state Potts 
models. Clearly the data of distinct lattice sizes collapse 
in a straight line in agreement with (j3)). The angular 
coefficient of these lines gives us the estimate c(L), for 
lattice size L. In table 1 we give these estimates for the 

TABLE I: Numerical values of the constant c(L) for Q — 
2, 3, 4 Potts models, the XXZ model (XXZ A ) and the Fateev- 
Zamolodchikov (FZ^) model. The expected values for the 
conformal anomalies together with the lattice sizes used in 
the numerical calculation are also shown. 
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FIG. 3: (Color online) Mutual information as a function of 
the subsystem size for the XXZ model with different values 
of the anisotropy A and lattice sizes L. 



expected value for the conformal anomaly for the model. 
The agreement with the predicted values are remarkable 
already for the lattice sizes we considered. 

XXZ quantum chain - The model describes the dy- 
namics of spin-i particles given by the Hamiltonian 

L 

Hxxz - - EK CT i+i + + A<7 K+i)' ( 9 ) 

1=1 

where a x ,a v cr z are spin-i Pauli matrices and A an 
anisotropy. This model provides an interesting check for 
the universal behavior of the Shanon mutual information. 
It has a continuous critical line, — 1 < A < 1, whose CFT 
has a central charge c = 1. According to ([3]) we should 
expect a data collapse of I(£,L) — I(L/2,L) for distinct 
lattice sizes and anisotropics. 

In order to illustrate this fact we plot in Fig. 3 the 
ratio [J(^,L)-/(L/2,L)]/[ln(sin(7r^/L))/4] as a function 
of £/L. We clearly see that for distinct anisotropics and 
lattice sizes the ratio is close [2l[ to a constant value 
given by the conformal anomaly. Similarly as we did 
for the Potts models the estimated values c(L) for the 
anisotropics A = and A = ±| are shown in table 1. 

Fateev-Zamolodchikov model - This is a spin 1 
model whose Hamiltonian is given by [22j 

L 

H FZ = - [of f - 2(cos 7 - l)(tr,V + of 

i=l 

-2sin 2 7 (af-(^) 2 +2(S 4 2 )2 }; (10) 

where S = (S x , S y , S z ) are spin-1 SU(2) matrices, <7j = 
SfS* +1 and cr.; = Sj-Sj+i = af- + erf. The model is anti- 
ferromagnetic for 6 = +1 and ferromagnetic for e = — 1. 
This is an important check for the conjecture © since 
the model has a line of critical points (0 < 7 < -|) with a 
quite distinct behavior in the antifcrromagnctic (e = +1) 
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FIG. 4: (Color online) Mutual information dependence with 
the subsystem size for the Fateev-Zamolodchikov model with 
different values of 7 and lattice sizes. 
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FIG. 5: (Color online) The difference I{£, L) - I(L/2, L), as a 
function of ln(sin(7r^/L))/4, for some lattice sizes for the XXZ 
chain with anisotropics A = —1/2 and A = 0, the ferromag- 
netic spin-1 Fateev-Zamolodchikov model with anisotropies 
7 = 7r/3 and 7 = 7r/4 and for the 4-state Potts model. All 
these models are ruled by a CFT with conformal anomaly 
c= 1. 



and ferromagnetic (e = — 1) cases. The antiferromagnetic 
version of the model is governed by a CFT with central 
charge c = | 12 3| while the ferromagnetic one is ruled by 
a c= 1 CFT HI. 

In Fig. 4 we show for e = ±1, the ratio [I(£,L) — 
/(L/2,L)]/[ln(sin(7r£/L))/4] as a function of tjL. The 
data are shown for the anisotropics 7 = 7r/3 and 7 = 7r/4 
and for lattice sizes L = 18 and L = 20. We clearly see 
an agreement with the expected central charge of the cor- 
responding CFT. In table 1 we give the predicted values 
c(L) obtained from the largest lattice size we considered 
for this model, and the agreement is again quite good. 

Conclusions - All the analytical and numerical cal- 



culations presented in this Letter indicate the following 
properties for the Shannon mutual information for the 
ground state of critical quantum chains: a) The leading 
dependence with the subsystem size I is independent of 
the basis in which we write down the ground state wave- 
function and characterizes the universality class of the 
critical behavior of the quantum chain, b) The finite-size 
scaling function ln(— sin ^) is the same as that of the en- 
tanglement entropy, c) The finite-size scaling, similarly to 
the entanglement entropy, is proportional to the central 
charge of the underlying CFT. An overall illustration of 
these points is presented in Fig. 5, where we show for 
2 < £ < L — 1 the finite-size behavior of the Shannon 
mutual information for the different models with central 
charge c = 1 presented in this letter. This figure shows 
that models whose Hamiltonians acts on rather distinct 
Hilbcrt spaces share the same universal behavior for the 
Shannon mutual information of their ground states. The 
numerical estimates presented in this Letter indicate that 
the Shannon mutual information as well as the entangle- 
ment entropy provide excellent tools for the evaluation of 
the central charge of conformal invariant quantum chains. 
It is important to mention that although the numerical 
results presented in this Letter are obtained by using the 
Lanczos method, we also verified that the Shannon mu- 
tual information can also be computed [25^ , for relatively 
large lattices, by using the DMRG (26j . 

Finally we mention that the conjecture (|3|) announced 
in this paper raises several interesting questions to be 
answered in the future. The first one concerns its proof 
based on CFT calculations as done by Calabresc and 
Cardy Q in the case of the entanglement entropy. In 
the case of the harmonic oscillator chain we show that 
in the bulk limit the dominant part of the Shannon mu- 
tual information is the same as that of the n = 2 Renyi 
entanglement entropy. A proof of this equivalence would 
produce as a corollary the conjecture (|3]). The Shannon 
mutual information of quantum chains with other bound- 
ary conditions as well as the study of the Rcnyi mutual 
information for the quantum chains are also interesting 
points to be studied in the future. The entanglement en- 
tropy of pure states formed by excited states also shows 
universal properties with finite size scaling functions re- 
lated to correlations of the associated CFtJBJ, the ques- 
tion whether this also happens for the Shannon mutual 
information is also an interesting one. As we can notice in 
Figs. 3 and 4, similarly as the entanglement entropy the 
Shannon mutual information also shows even-odd parity 
oscillations with the subsystem size I J^- The study 
of these oscillations is a also an interesting point for the 
future. 

This work was supported in part by FAPESP and 
CNPq (Brazilian agencies). We thanks J. A. Hoyos and 
R. Pereira for useful discussions and a careful reading of 
the manuscript. 



5 



[1] L. Amico,R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. 

Phys. 80, 517 (2008). 
[2] K. Modi, A. Brodutch, H. Cable, T. Paterek, and V. 

Vedral, Rev. Mod. Phys. 84, 1655 (2012). 
[3] C. Holzhey, F. Larsen, and F. Wilczek, Nucl. Phys. B 

424, 443 (1994). 
[4] P. Calabrese and J. Cardy, J. Stat. Mech. (2004) P06002 

[5] F. C. Alcaraz, M. I. Berganza, and G. Sierra, Phys. Rev. 
Lett. 106, 201601 (2011); M. I. Berganza, F. C. Alcaraz, 
and G. Sierra, J. Stat. Mech. (2012) P01016; L. Taddia, 
J. C. Xavier, F. C. Alcaraz, and G. Sierra (unpublished). 

[6] J. Cardy, Phys. Rev. Lett. 106, 150404 (2011). 

[7] D. A. Abanin and E. Demler, Phys. Rev. Lett. 109, 
020504 (2012). 

[8] A. J. Daley, H. Pichler, J. Schachenmayer, and P. Zoller, 

Phys. Rev. Lett. 109,020505 (2012). 
[9] J. Urn, H. Park and H. Hinrichsen, J. Stat. Mech. (2012) 

P10026 . 

[10] M. M. Wolf, F. Verstraete, M. B. Hastings, J. I. Cirac, 

Phys. Rev. Lett. 100, 070502 (2008). 
[11] H. Bernigau, M. J. Kastoryano, and J. Eisert, 

[arXiv: 1301.5646| . 
[12] J-M Stephan, S. Furukawa, G. Misguich, and V. 

Pasquier, Phys. Rev. B, 80, 184421 (2009). 
[13] J-M Stephan, G. Misguich, and V. Pasquier, Phys. Rev. 

B, 82, 125455 (2010); Phys. Rev. B 84, 195128 (2011). 
[14] M. P. Zaletel, J. H. Bardarson, and J. E. Moore, Phys. 

Rev. Lett. 107, 020402 (2011). 
[15] J. Wilms, M. Troyer, and F. Verstraete, J. Stat. Mech. 



(2011) P10011; J. Wilms, J. Vidal, F. Verstraete,and S. 
Dusuel, J. Stat. Mech. (2012)P01023. 
[16] H. W. Lau and P. Grassberger, Phys. Rev. E 87, 022128 
(2013). 

[17] R. G. Unanyan and M. Fleischhauer, Phys. Rev. Lett. 

95, 260604 (2005). 
[18] M. Cramer, J. Eisert, M.B. Plenio and J. Dreifiig, Phys. 

Rev. A 73, 012309 (2006). 
[19] H. Casini and M. Huerta, J. Phys. A 42, 504007 (2009). 
[20] F. Y. Wu, Rev. Mod. Phys, 54, 235 (1982). 
[21] The small deviations that appears for £ = L/2 — 1 is 

enlarged due to the division by ln(sin(7r£/L)). 
[22] A. B. Zamolodchikov and V. A. Fateev, Yad. Fiz. 32, 581 

(1980) [Sov. J. Nucl. Phys. 32, 298 (1980)]. 
[23] F. C. Alcaraz and M. J. Martins, Phys. Rev. Lett. 61, 

1529 (1988); P. Di Francesco, H. Saleur, and J.-B. Zuber, 

Nucl. Phys. B 300, 393 (1988). 
[24] F. C. Alcaraz and M. J. Martins, Phys. Rev. Lett. 63, 

708 (1989). 

[25] F. C. Alcaraz, M. Rajabpour, and J. C. Xavier (unpub- 
lished) . 

[26] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

[27] P. Calabrese, M. Campostrini, F. Essler, and B. Nienhuis, 
Phys. Rev. Lett. 104, 095701 (2010); P. Calabrese and F. 
H. L. Essler, J. Stat. Mech. (2011)P08029; M. Fagotti and 
P. Calabrese, J. Stat. Mech. (2011)P01017; J. Cardy and 
P. Calabrese, J. Stat. Mech. (2010)P04023; J. C. Xavier 
and F. C. Alcaraz, Phys. Rev. B 85, 024418 (2012). 



